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ON THE GIBBS PHENOMENON I: 

RECOVERING EXPONENTIAL ACCURACY FROM THE 
FOURIER PARTIAL SUM OF A NON-PERIODIC ANALYTIC FUNCTION 1 


David Gottlieb, Chi-Wang Shu, Alex Solomonoff and Herve Vandeven 2 
Division of Applied Mathematics 
Brown University 
Providence, RI 02912 

ABSTRACT 

It is well known that the Fourier series of an analytic and periodic function, truncated 
after 2 N + 1 terms, converges exponentially with V, even in the maximum norm. It is 
also known that if the function is not periodic, the rate of convergence deteriorates; in 
particular there is no convergence in the maximum norm, although the function is still 
analytic. This is known as the Gibbs phenomenon. In this paper we show that the 
first 2N + 1 Fourier coefficients contain enough information about the function so that an 
exponentially convergent approximation (in the maximum norm) can be constructed. The 
proof is a constructive one and makes use of the Gegenbauer polynomials C^(x). It consists 
of the following two steps. 

In the first step we show that the first m coefficients of the Gegenbauer expansion (based 
on C^(x), for 0 < n < m) of any L 2 function can be obtained, within exponential accuracy, 
provided that both A and m are proportional to (but smaller than) N . 

In the second step we construct the Gegenbauer expansion based on C* , 0 < n < m 
from the coefficients found in the first step. We show that this series converges exponentially 
with N, provided that the original function is analytic (thou gh non-periodic). 

Thus we prove that : The Gibbs phenomenon can be completely overcome. 
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1 Introduction 


In this paper, we deal with a prototype of the Gibbs phenomenon and show how to eliminate 
it. Consider an analytic but non-periodic function f(x) defined in [—1,1]. Notice that f(x) 
has a discontinuity at the boundary x = ±1 if it is extended periodically with period 2. The 
Fourier coefficients of /(x) are defined by 

/(*) = \ J\ . ( 1 . 1 ) 

Assume that the first 2N + 1 Fourier coefficients /(fc), |fc| < N, are known but the 
function /(x) is not. Our objective is to recover the function /(x) for —1 < x < 1 with 
exponential accuracy in the maximum norm. The traditional Fourier partial sum using the 
first 2N + 1 modes 


!n(x)= £ /(*)«“” (1.2) 

k=-N 

does a poor job; it produces a first order approximation to f(x) with an error O(^) away 
from the boundary x = ±1, and shows 0(1) spurious oscillations near the boundary x = ±1 
known as the Gibbs phenomenon. Thus there is no convergence in the maximum norm. 
When one uses a filter in the Fourier space, 

/»(*) = E of /Me*” , (1.3) 

k=-N 

where a £ = a^ k are suitably defined real or complex numbers which tend to zero when \k\ 
tends to N , the situation becomes better. One can get exponential accuracy away from the 
boundary x = ±1 if a £ are chosen as suitable real numbers [7, 9, 10, 12] or one can get 
exponential accuracy up to one boundary x = — 1 or x = 1 if aj? are chosen as suitable 
complex numbers [3]. In these cases the approximation is still in the space spanned by the 
first 2N + 1 trigonometric polynomials and is a convolution of the original Fourier partial 
sum with some filter kernel which is an approximate two-sided or one-sided S function; hence 
it cannot be exponentially accurate in the maximum norm for — 1 < x < 1. For the one- 
sided filters introduced in [3] one can use two different approximations in — 1 < x < 0 and 
in 0 < x < 1, right-sided for the former and left-sided for the latter, to obtain exponential 
convergence globally. 

In this paper we adopt a different point of view. The idea is the following: we realize 
that the problem with the Fourier approximation is the non-periodicity of the function and 
the fact that the functions e mkx are the solutions of a regular Sturm-Liouville problem. In 
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[6] it is shown that expanding an analytic, non-periodic function f(x) by the eigenfunctions 
of a singular Sturm-Liouville problem yields rapid convergence. For example, a Chebyshev 
or Legendre expansion of f(x) converges exponentially. Thus if the first 2N + 1 Fourier 
coefficients can provide enough information to reconstruct the coefficients of an expansion 
based on a singular S-L problem, we might recover the accuracy. Unfortunately, one can 
not recover the coefficients of the Chebyshev or the Legendre expansion within high enough 
accuracy. 

In this paper we show that from the first 2iV + 1 Fourier coefficients of an analytic but 
non-periodic function, one can get the first m ~ N coefficients in the Gegenbauer expansion 
based on the Gegenbauer polynomials C„(x), provided that the parameter A, appearing in 
the weight function (1 — x 2 ) x ~i , grows with the number of Fourier modes- N . We prove that 
this yields exponential accuracy in the maximum norm. 

Our proof consists of two separate and independent steps. The first step (Section 3) is to 
show that given the Fourier partial sum of the first 2N+1 Fourier modes, of an arbitrary L 2 
function f{x), it is possible to recover the partial sum of the first m terms in the Gegenbauer 
expansion of the same function to exponential accuracy (in the maximum norm) by letting 
the parameter A and the number of terms m in the Gegenbauer expansion grow linearly with 
N. In this step f(x) need not be smooth; any L 2 function will do. We denote the error 
between the exact Gegenbauer coefficients and those obtained from the Fourier coefficients 
- the Truncation error. The results of this Section are summarized i n Theorem 3.3 . 

In the second step (called the Regularization error), we prove the exponential convergence, 
in the maximum norm of the Gegenbauer expansion of an analytic function when A grows 
linearly with m. This is done in Section 4. The second step has its own interest; it is 
an exponential convergence proof in the maximum norm for such Gegenbauer expansions of 
analytic functions, where A increases with the number of the terms used in the approximation. 
The results of this Section are summarized in Theorem 4.3 . 

In Section 2 we discuss some results concerning Gegenbauer polynomials that are relevant 
to the proof and the computations. These are combined with the results of Sections 3 and 
4 to yield the main theorem of Section 5 which demonstrates that one can construct an 
exponentially convergent approximation to an analytic, non-periodic function, from its first 
2N + 1 Fourier coefficients. Finally, in Section 6 we demonstrate the theory with some 
numerical examples. Of special interest is Example 6.1, concerning the function f(x) = x. 
This function was used originally (in 1898) to demonstrate the Gibbs phenomenon. 

We will use A or A for a generic constant independent of all the growing parameters 
throughout this paper. The actual value of A or A may be different in different locations. 
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2 Preliminaries 

In this section we will introduce the Gegenbauer polynomials and discuss some of their 
asymptotic behavior. We rely heavily on the standardization in Bateman [2], although for 
our purpose a different scaling might have been more natural. 

We start by defining the Gegenbauer polynomials C*(x) in the following definition. 

Definition 2.1 

The Gegenbauer polynomial C*(x) is the polynomial of order n that satisfies 


J^(l-x 2 ) x -$C£(x)C!;{x)dx = 0 k^n 


and (for A > 0) 


A _ F ( n + 2A ) . 

n!r(2A) 


( 2 . 1 ) 


(2.2) 


Note that the Gegenbauer polynomials thus defined are not orthonormal. In fact, the 
norm of C*(x) is given by the following lemma. 


Lemma 2.1 

The Gegenbauer polynomials defined above satisfy 


j\( 1 - x 2 ) x -iC^(x)C^x)dx = hi 

where 



ci(i) 


r(A + 

r(A)(n + A) ' 


For the proof see [2, page 174]. 


(2.3) 


(2.4) 


□ 


We are ready now to deal with the asymptotics of the Gegenbauer polynomials for large 
n and A. For this we need Lemma 2.2 and 2.3. 

Lemma 2.2 (Stirling) 

For any number x such that x > 1 we have 

(2n)h x+ h- x eh 

{2ir)h x+ h' x 


r(* + i) < 

r(x + i) > 


3 


(2.5) 

( 2 . 6 ) 


Lemma 2.3 

There exists a constant A independent of A and n such that 


□ 


The proof follows from (2.4) and the Stirling’s formula (2.5)-(2.6). 


(2.7) 

( 2 . 8 ) 


□ 


Finally we would like to quote the Rodrigues ’ formula [2, page 175]. 


Lemma 2.4 

The Gegenbauer polynomials are explicitly given by 


(1 - = U £ g ( A ,„)^ [(1 - 

where G(\,n ) is defined by 


G{X,n) 


T(A + i)r(n + 2A) 
r(2A)r(n + A + 1) • 


(2.9) 


( 2 . 10 ) 

□ 


3 The Truncation Error 

In this section we consider an arbitrary L 2 function f(x ) defined in [—1,1]. We assume that 
the first 2 N + 1 Fourier coefficients f(k), as defined in (1.1), are given. We are interested 
in recovering (within exponential accuracy in the maximum norm) the first m coefficients in 
the Gegenbauer expansion of f(x), 

OO 

/(*) = £ /*(')<?(*), (3.1) 

l - 0 

where the Gegenbauer coefficients are defined by 
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(3.2) 


/'*(/) = Jj jj 1 - x ! ) x ^/(x)C, x (x)dx 

with h x given by (2.4). 

Since we do not know the function f(x), but rather its truncated Fourier series /n(x) 
(defined in (1.2)), we have only an approximation to f x (l) which we denote by g x (l), given 
by 


S x (l) = Jj 1 - xV-t Mx)C,\x)dx . 


(3.3) 


Notice that g x (l) depends on N. At this stage we would like to define the truj^cgjjon error 


TE(X,m,N) 


max 

-l<r<l 


E(/ A (0- » A (0)c A M . 


/=0 


(3.4) 


The truncation error is the difference between the Gegenbauer expansion (with m terms) 
of the function f(x) and that of the truncated Fourier series /n(x), It measures the error in 
the finite Gegenbauer expansion due to truncating the Fourier series. 

In the next two theorems we bound the truncation error in terms of N, the number of 
given Fourier coefficients, m the number of Gegenbauer polynomials, and A. 


Theorem 3.1 

If f(x) is an L 2 function on [—1, 1], then there exists a constant A which is independent 
of A, m and N, such that the truncation error defined in (3.4) satisfies the following estimate: 


TE(\,m,N ) 


^ (m + A)T(m + 2A)r(A) 
(m — l)!r(2A) 



(3.5) 


Proof: 

As a first step we consider the special function f(x) = e ,nnx with |n| > N. In this special 
case /n(x) = 0 and we obtain 

(/*(!) - j A (0)C, J (l) = jj 1 - x J ) A -l e '“"C, A (x)^ . (3.6) 

Roughly speaking, one can argue that this integral is rapidly decreasing when n increases. 
This results from the fact that the integral is proportional to the n’th Fourier coefficient of 
the function (1 — x 2 ) x ~^C x (x) which is analytic and has A periodic derivatives. It is nice to 
know that an explicit expression of this integral appears in the literature [2, page 213]: 
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(3.7) 


J (1 - x 2 ) A *e ,nwx Cf(x)dx = T(A) i‘(l + A)J; + a(ttw) 

where J„(x) is the Bessel function. Since |J^(a:)| < 1 for all x and v > 0 [1, page 362], we 
have, for 0 < / < m, 

- g\l))CKl)\ < 


< 

where in the second step we used the formula (2.2) for C A (1), and in the last step we used 
the fact that ( , +*)H i + 2 ' x ) j s an increasing function of /. 

We now return to the general function /(x), which satisfies 


c '<w a >Uk) (,+A) 

(M-A)r(l + 2A)r(A) / 2 A a 
/! r(2A) \ 7 r|n|y 

(m + A)r(m + 2A)r(A) / 2 V 
m!r(2A) \ 7r l n l/ 


(3.8) 


f(x) - f N (x) = E /(«) e * W • 

|n|>AT 

Since f(x) is an L 2 function, its Fourier coefficients f(n) are uniformly bounded, 


(3.9) 


l/(«)l < A ■ 

We thus have, using the result for the special case e' nirx in (3.8), 

(m + A)r(m + 2A)r(A) 


(3.10) 


l(/ A (0- $*(0)<?(1)I < A 

< A 


m\ T(2A) 

(m + A)r(m + 2A)r(A) / 2 ' A " J 
m!r(2A) 


£(i) 


(3.11) 


for all 0 < / < m. 

We can now estimate the truncation error (3.4) by 


TE(\,m,N ) < 

< 

< 


m max max |(/ A (/) - £ A (/))C A (x)| 

0<l<m - 1 < x <1 

m max |(/ A (0 - 9 X (l))\C?{l) 

0 </<m 


<71 A 


(m + A)r(m + 2A)r(A) ( 2 \ A_1 


m 


ir(2A) 


\xN) 


(m + A)r(m + 2A)r(A) / 2 n"' 1 
(m - l)!r(2A) \irNj 


6 



where in the second step we used the fact that IC^z)! < C;( 1) for all — 1 < x < 1 
[2, page 206], and in the third step we used (3.11). 

The theorem is now proven. 


□ 


Remark 3.1 The approximate Gegenbauer coefficients g > '( l ) in (3.3) can be explicitly 
expressed in terms of the Fourier coefficients f(k ) as 

g\l) = W(0) + r(A)i'(/ + A) £ 4 +a(»*)( 4 ) A /(*)- (3.12) 

0<|fc|<7V \KK/ 

Equation (3.12) follows immediately from the definition (3.3) of g*(l) and the integration 
formula (3.7). 


□ 

For fixed A, the truncation error (3.5) decays algebraically as O ( N x-i )• However, if both 
A and m grow linearly with N, the truncation error can be made exponentially small. This 
is stated in the following theorem. 

Theorem 3,2 

If A = aN and m = (3N where a and 0 are positive constants, then the truncation error 
defined in (3.4) satisfies 


TE(aN, 0N, N ) < AN 2 q N (3.13) 

where 


(g + 2 

q (2tt c)"q“/J 0 ’ <3 l4) 

In particular, if a = then 


q = e 27 « 0.8 < 1 . 


Proof; 

We use Stirling’s formula (2.5)-(2.6) to obtain, from the previous estimates on the trun- 
cation error in (3.5) and some simple algebra, 
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TE(aN,/3N,N) 


< 

< 


A{/3 + a)N 


T((fi + 2a)N)T(aN) 
1 0N - l)\T(2aN) 


AN 2 q N 


/ 2 \ °>N-\ 

W/vJ 


with q defined by (3.14). If we take a = j3 in (3.14), we obtain 



which attains its minimum value q = e ™ at /? = ||. 


□ 


We would like to point out that we choose a = j3 in Theorem 3.2 simply to show that it 
is possible to obtain exponentially small truncation errors. This may not be the best choice 
in practice. We can easily verify that, for fixed a, (3.14) defines a q which is an increasing 
function of /?. This is not surprising since the truncation error should be bigger if there are 
more terms in the Gegenbauer expansion to approximate. However we will see in the next 
section that the regularization error will be smaller if m is bigger. In practice one might try 
to choose m to maintain some balance between these two errors. For a fixed (3 , (3.14) defines 
a q which attains its minimum at 

tt — 2/3 + Jt^k - 4/3) 

a = p (3.15) 

if f - For example, if /? = |y and a « 1.33 is chosen according to (3.15), then q given by 
(3.14) is approximately 0.49, much smaller than the minimum value 0.8 obtained with the 
restriction a = j3. 

We summarize the results of this Section in the next theorem. 

Theorem 3.3 The exponential decay of the truncation error 

Let f(x) be an L 2 [— 1, 1] function, and f(k), —N < k < N, its Fourier coefficients defined 
in (1.1). Let f x (l) be the Gegenbauer expansion coefficients of f(x ) defined in (3.2), and let 
g x (l ) be the Gegenbauer coefficients of the truncated Fourier series /n(x) defined in (1.2) 
with g x {l) are given explicitly in (3.12). 

Then: if A = m = (3N, where j3 < ^r, the truncation error decays exponentially with 
the number of Fourier modes N , i.e. 

m 

TE(/3N , fiN t N ) = max | £(/ A (/) - g x (l))C x (x)\ < AN\ n (3.16) 

1=0 
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with 


9 


( 27£) f 

^27re' 


< 1 


□ 


4 The Regularization Error 

In this section we would like to establish error estimates for approximating an analytic 
function f(x) on [—1, 1] by its Gegenbauer expansion based on the Gegenbauer polynomials 
C£(x). Since our goal is to remove the Gibbs phenomenon we will use the maximum norm. 
In the last section we have shown that we can get the Gegenbauer partial sum of the first 
m terms of any L 2 function from its Fourier partial sum of the first 2 N + 1 modes with 
exponential accuracy in the maximum norm, if A and vn are both growing linearly with N . 
Thus in this section we will consider the case of large A and m. 

We will assume that f{x) is an analytic function on [-1,1] satisfying the following as- 
sumption. 


Assumption 4.1 

There exists constants p > 1 and C(p) such that, for every k > 0, 


max 

-l<r<l 





(4.1) 


This is a standard assumption for analytic functions. Here p is actually the distance from 
[-1,1] to the nearest singularity of f(x ) in the complex plane (see for example [8]). The 
assumption can be modified using the techniques in [5]. 

Let us consider the Gegenbauer partial sum of the first m terms for the function f(x ) 
given by (3.1), with the Gegenbauer coefficients f x {l) defined by (3.2). We want to estimate 
the regularization error in the maximum norm, 


RE(\,m).= max 

V ' -Kx< 1 


rn 


1 = 0 


We start by estimating the Gegenbauer coefficients f x (l)- 


(4.2) 


Lemma 4.1 

The Gegenbauer coefficient / A (/), as defined in (3.2), of an analytic function satisfying 
the Assumption (4.1), can be bounded by 
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lf(OI < A 


CM r(A + OT + 2A) 
A, A (2/>)T(2A)r(/ + A + 1 ) ' 


(4.3) 


Proof; 


We start by using the definition (3.2) for / A (/). We replace the term (1 - x 2 ) x ~i C x (x) 
by the Rodrigues’ formula (2.9)-(2.10) to get 

m) = /<«>sr [d - * 

where G{ A, 1) is defined in (2.10). Integrating by parts l times we get 

We now use the Assumption (4.1) to estimate the derivative ^f(^), thus obtaining 

|/ A (0I < /‘(l - x 2 )' +A -trfx . 

Since Cq (x) = 1, the remaining integral is simply /io +A and can be obtained from (2.3)-(2.4): 

I fx,n\ < G(A,/)C(p)V?T(/ + A + 1) 

l/U| - ^ 2 V(/ + A)r(/ + A) 

and finally using the definition of G(\,I) from (2.10) we get (4.3). 


The estimate (4.3) can be used naturally to get an estimate in the weighted L 2 norm. 
This would have been more transparent if we had adopted an orthonormal Gegenbauer basis 
rather than (2.2). However the aim of this paper is to establish estimates in the maximum 
norm. This is stated in Theorem 4.1. 


Theorem 4.1 

If f(x) is an analytic function on [-1, 1] satisfying the Assumption (4.1), then the regu- 
larization error defined in (4.2) can be bounded by 


RE(X,m) < ^g O’)r(* + j)r(™ + 2A+l) 

“ m\/A(2p) m r(2A)r(m + A) 

Proof: 

By (4.3) and (2.8) we obtain 


(4.4) 
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(4.5) 


If we define 


l/ A (Oltf(l) < A 


c(<.)r(A + |)r(i + 2 A) 
\/A(2/j) , r(2A)r(/ + A) ' 


nm _ , cwr(A + i )r (H-2A) 

' ’ v / A(2/))'r(2A)r(l + A) 


then clearly for m < / 

B(l+ 1) 1 + 2A l + £ < 1 _+jt 

5(/) 2/,(/ + A)-p(2 + f)-2 + ^' 


We can thus sum (4.5) for m + l</<ooto obtain 


RE(\,m ) = 


max ^ f x (l)C?(x) 


i=m+1 


< E l/ A (')l<7(i) 

/=m+l 

oo 

< E *(') 

/=m-f 1 

,2( m + A) 

m 

C(p) r(A + *)r(m + 2A + 1) 
m\/A(2/?) m r(2A)r(m + A) 


< 5(m + l)- 

< A- 


where we have used again the fact that |C*(x)| < Cf ( 1 ) for all 1 < x < 1 in the second 
step. This finishes the proof. 


□ 


About the size of the regularization error (4.4) when A depends linearly on m, as is the 
case in Theorem 3.2 for truncation errors, we can state the following. 


Theorem 4.2 

If A = 7m where 7 is a positive constant, then the regularization error defined in (4.2) 
satisfies 


55(7m,m) < Aq m 


(4.6) 


where q is given by 
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which is always less than 1. 
constant, then 


= (1 + 27) x +^ 
q /?2 1 + 2 t' 7T'(1 +7) 1 +'r 
In particular, if 7 = 1 and m = 


(4.7) 


/3N where is a positive 


RE(PNJN) < Aq N 

with 


(4.8) 



Proof: 


(4.9) 


We use Stirling’s formula (2.5)-(2.6) to replace the Gamma functions in (4.4). A little bit 
of algebra brings us to (4.6)-(4.7) if A = 7 m. Notice that the constant A in (4.6) contains 
the contribution of p related terms. It is easy to verify that q defined by (4.7) is an strictly 
increasing function of 7 and tends to A < 1 when 7 tends to infinity. Hence we have q < A < 1 
for all 7 > 0. As for the proof of (4.8)-(4.9), we simply plug in 7 = 1 and m = /3N into 
(4.6)-(4.7). 


□ 


Finally we summarize the results of this Section in the following theorem. 

Theorem 4.3 The exponential decay of the regularization erro r 

Let f(x) be an analytic function on [—1, 1] satisfying assumption (4.1). Let 0 < / < 
m be its Gegenbauer coefficients defined in (3.2). Furthermore assume that A = m = 0N. 
Then 


max 

-1<X<1 


/(*) - E /*(')#(*) 


l=o 


<Aq N 


(4.10) 


where 



□ 
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5 The Main Theorem 


In this Section we bring the main theorem demonstrating that one can construct an ex- 
ponentially convergent (in the maximum norm ) approximation to an analytic, non-periodic 
function, from its first 2 N -1- 1 Fourier coefficients. The method is indicated in the last two 
sections. Namely, first we get, from the Fourier coefficients, an approximation to the first 
m - pN Gegenbauer coefficients. This approximation, by virtue of the discussion in Sec- 
tion 3, converges exponentially f ast to the true coefficients, provided that A grows with N . 
We then construct the partial Gegenbauer expansion which converges exponentially to f(x) 
by virtue of the discussion in Section 4. 

We are ready to state our main theorem. 

Theorem 5.1 Removal of the Gibbs Phe nomenon 

Consider an analytic and non-periodic function f(x) on [—1,1], satisfying 


max 

-1<X<1 





P > 1 • 


(5.1) 


Assume that the Fourier coefficients 

/t* o = \ £ 


are known for —N < k < N. 

Let g x {l), 0 < / < m be the Gegenbauer expansion coefficients of Jn(x) = J2k=-N f(k)e' k * x 
explicitly given by 


g\l ) = 6oif( 0) + r(A)i'(/ + A) £ ( 4 :) /(*) • ( 5 - 2 ) 

0<|Jt|<JV X7I "' / 


Then for A = m = 0N where 0 < we have 


Here 


max 

-1<I<1 


f( x ) ~Y,9 X ( 0C| A (*) 


l - 0 


< AN 2 qx + Mr 


q T = 




(5.3) 


Proof: 

We start by noting that formula (5.2), for the approximate Gegenbauer coefficients g x (l), 
follows from Remark 3.1. 

In order to establish (5.3), we introduce the Gegenbauer coefficients of the function f{x) 
and denote them by f x (l), they are defined in (3.2). 
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We have 


max 

-1<X<1 


m 


/=0 


< max 

“ -l<x< 1 


max 

— 1<x<1 


/(*)-£/*«<?(*) 


1—0 


+ 


£/*(<)<?(*) -E^rocfw 

i=0 (=0 


The first term is the regularization error and has been estimated in Theorem 4.3. The 
second term is the truncation error and has been estimated in Theorem 3.3. The theorem is 
thus proved. 


□ 


A few remarks are in order: 

Remark 5.1 The proof is constructive. Given 2 N + 1 fourier coefficients f(k), 
—N < k < N, one gets explicitly the approximate Gegenbauer coefficients g x (l) , 0 < l < m 
and the Gegenbauer series can be explicitly constructed. 

Remark 5.2 The reconstruction is not optimal, since no effort has been spent to optimize 
the parameters. 

6 Numerical Results 

In this section we demonstrate the theory using numerical examples. We implement the 
method in the following way. Assume that the first 2 N + 1 Fourier coefficients of f(x), f(k) 
for — N < k < N, as defined in (1.1), are given. We compute the approximate Gegenbauer 
coefficients g x (l), for 0 < / < m, defined in (3.3), using the following formula given in Remark 
3.1: 


= W(0) + r(A)i'(f + A) £ Ji+a(it*)(4) A /(*). (S- 1 ) 

0<|Jt|<JV 

We compute the Bessel function J„(x) using an IMSL routine. Once the approximate 
Gegenbauer coefficients g x (l) are obtained, we can compute the approximation to f(x) by 
directly summing 


»£(*) = E$ l (0C/(x) (6.2) 

/=0 

as long as we can compute the Gegenbauer polynomial C x (x) accurately. The formula we 
used to compute C x (x) is 
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(6.3) 


Cl' (cos 8) 


i 


E 


T(k + A) T(l-k + X) 
jbir(A) (/-Jt)ir(A) 


cos(/ — 2 k)8 


which can be found in [1, page 175]. 

We remark that the implementation described above are subject to round off effects for 
large A and m. We use a Cray-YMP to carry out all the computations. Our implementation 
can give accurate result only when the error is no smaller than 10 9 . And we will show results 
only in those cases. A better way to implement this method might be through Chebyshev 
polynomials. 

For the purpose of testing we consider the following examples. 


Example 6.1 f(x) = x 

This is the original example used in 1898 to show the famous Gibbs phenomenon. For 
this particular function there is no regularization error as long as m > 1, hence all the errors 
observed result from the truncation error. 

In Figure 1 , we show the exact function /(x) — x in the solid line, the Fourier partial sum 
f N (x) with N = 4 in the short dashed line, and the approximation through the Gegenbauer 
polynomials g^(x), as defined in (6.2), with N = 4 and m = A = \N in the long dashed 
line. We can clearly see that the Fourier partial sum f N (x ) shows the Gibbs oscillations, 
while the approximation gl n (x) through the Gegenbauer polynomials is uniformly accurate. 
In Figure 2 we draw the errors in the maximum norm, with a logarithmic scale, versus N 
(recall that 2N + 1 is the number of Fourier modes given), both with m — A = \N (squares) 
and with m = \N and A determined by (3.15) (circles). This picture confirms our estimates 
in (3.13) through (3.15) for the exponential convergence of the truncation errors. 

Next, we show the point-wise errors of the approximations with m = A = \N for 
AT = 4, 8, 16,32 in Figure 3, and the point- wise errors with m = \N and A determined 
by (3.15), for N = 4, 8, 16, in Figure 4, again both in logarithmic scales. We can observe ex- 
ponential convergence both inside the interval and near the boundary, although the absolute 
error is several magnitudes smaller inside the interval than at the boundary. 

We have also run the test for f(x) = x n with n = 2 through 9, obtaining similar results 
(not shown here). 

Example 6.2 f(x) = cos[1.47r(x - 1)]. 

This function satisfies the Assumption (4.1) with arbitrary p > 1. Both truncation 
error and regularization error will be present. The choice m = A = \N gives exponential 
convergence as expected. The choice m = A = \N gives better results in this case. 
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In Figure 5 we draw the errors in the maximum norm versus N , with a logarithmic scale, 
both with m = A — | N (squares) and with m — A = (circles). The choice of m = ^ N 
and A determined by (3.15) does not converge in this case, indicating that (3.15) is probably 
a bad choice for the Regularization errors. 

Finally, again in logarithmic scales, we show the point-wise errors of the approximations 
with m = A = | jV for N = 24, 36, 52, 76 in Figure 6, and the point-wise errors with m = 
A = | N for N = 16, 24, 36, 52 in Figure 7. We again observe exponential convergence both 
inside the interval and at the boundary, and the several magnitudes difference in the absolute 
errors inside the interval and at the boundary. 


7 Concluding Remarks 

We considered in this paper a prototype of the Gibbs phenomenon and demonstrated that 
this effect may be overcome. We show that if f(x) is an analytic but non-periodic function 
then the knowledge of its first 2 N + 1 Fourier coefficients provides accurate enough infor- 
mation so that the first m ~ N coefficients in the expansion of f(x ) by the Gegenbauer 
polynomials C£(x) (A ~ N ) may be obtained. In equation (5.2) the formula that relates the 
Gegenbauer coefficients to the Fourier coefficients is presented. 

There are several topics which are not addressed in this paper and will be discussed in a 
future paper: 

(1) Other choices of polynomials such as the Laguerre polynomials and Hermite polyno- 
mials; 

(2) Algebraic convergence for C k but not analytic functions; 

(3) The case with exponential recovery in a sub-interval [a, 6] of [—1, 1], That is, given 
the first 2 N + 1 Fourier modes for a function defined on [-1,1], find its Gegenbauer partial 
sum of the first m terms based on the scaled Gegenbauer polynomials on [a,b], in which 
the function f(x) is analytic. This would allow one to handle multiple discontinuities and 
discontinuities of unknown location; 

(4) The procedure realized in the Fourier space as a convolution, similar to the approach 
used in [3]; 

(5) Efficient and stable numerical implementation of those methods; 

(6) Extension to Collocation techniques and other than the Fourier basis functions. 
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True and Approximate Solutions 





Figure 2: For the function f(x) = x, the errors in the maximum norm, max-^o \f{x) - 
( x ) | where g^(x) is defined in (6.2), in a logarithmic scale, versus N{2N + 1 is the number 
of Fourier modes given). The squares are for the case with m = A = \N , and the circles are 
for the case with m = \N and A determined by (3.15). 







Logarithmic Errors in Maximum Norm 


Figure 5: For the function f(x ) = cos[1.47r(x — 1)], the errors in the maximum norm, 
max_i< 3: < 1 |/(x) — g^ (x) | where g^{x) is defined in (6.2), in a logarithmic scale, versus 
N (2 N + 1 is the number of Fourier modes given). The squares are for the case with 
m = A = ifV; and the circles are for the case with m = A = %N. 
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